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Abstract: Formulae are derived for the positions of end-points in the invariant mass 
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1. Introduction 

In recent years there has been considerable theoretical and experimental interest in the 
possibility of observing and measuring the production of heavy neutral SUSY higgs particles 
{H/A) at the LHC. In many regions of parameter space decays to Standard Model (SM) 
particles such as H/A —> tt provide the clean signature of a peak at m(H/A) in the 
reconstructed invariant mass distribution of the decay products. In certain regions of 
parameter space however, particularly at larger values of tra, this channel is suppressed and 
alternative signatures must be sought. One interesting possibility |l], |], ||, [|, |5|, [| |7|, ^, [l(J 



11, 12, [L3|] involves the use of higgs decays to pairs of gauginos (neutralinos or charginos), 
with subsequent gaugino decays to dileptons and the Lightest Super symmetric Particle 
(LSP - here assumed to be neutral and stable). The resulting 4£ + E™ 1SS signature suffers 
from relatively little irreducible Standard Model background contamination, although there 
can be significant irreducible SUSY background from direct or indirect gaugino production 
(the latter via squark decays). 

Several detailed studies have been performed of the likely sensitivity of the LHC de- 
tectors to SUSY higgs production using the 4£ + E™ 1SS signature (see e.g. Refs. Q [2], |3|, 
IE), ||, [?], ||, ||, |I0, |ll], [l2], 13 1). If a statistically significant excess is observed in this channel 



then it is natural to ask whether the same channel can be used to measure the mass of the 
higgs. Unlike in the case of SM decays it is likely to be extremely difficult to reconstruct a 
higgs mass peak due to the presence of the two heavy invisible LSPs in each higgs decay. 
In the case where the decay proceeds through the chain: 

h/a -> x° 2 x°2 -»• tM mixlxl (1.1) 
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and the masses of the X2, ^ an d Xi (LSP) are known, it is possible to solve the kinematic 
constraints to obtain an event-by-event value of m(H/A) ||. If the sleptons I are heavy 
however, causing the \\ particles decay through the three-body process X2 ~~ ^ ^Xi-> there 
are insufficient constraints to solve for m(H/A). In this case one must resort to measuring 
the positions of end-points in the distributions of event invariant mass and transverse mass 



|14, 15, [l^, 17, 18, 19, [2(| values, which depend on the masses of all the particles involved 
in the decay, including m(H/A). 

In this brief paper we shall derive formulae for the positions of these end-points both 
in the special case where the decay chains are identical and in the more general case of 
non-identical decay chains. The resulting expressions will be applicable to any channel in 
which a heavy state u decays to a pair of lighter states 5i (i = 1,2), each of which decays 
in turn to an aggregate visible state V{ and an invisible particle a{\ 

uj — > 8182 — > v\a\V20L2- (1.2) 

Note that the aggregate visible states Vi can each be composed of multiple visible particles 
- the V{ should be considered to be pseudo-particles with four-momenta equal to the net 
four-momenta of their constituents. We shall test the new formulae with a simple case- 
study of heavy SUSY higgs decays to gauginos. We shall not attempt to perform a full 
detector-level study of the SM and SUSY backgrounds to the AH, + E™ ms channel as this 
is described elsewhere (see e.g. Refs. [f|, |j| |l^, [OJ), however we shall discuss the relative 
merits of using the invariant mass and transverse mass end-points to measure m(H/A). 

The structure of the paper is as follows. In Section Q we derive the invariant mass end- 
point formulae while in Section [3| we derive the equivalent formulae for transverse mass 
end-points. In Section ||| we present the results of the case-study of SUSY higgs decays to 
four leptons and E^ 1SS . In Section || we conclude. 



2. Invariant mass end-points 

As is well-known the invariant mass of the aggregate visible products of each chain, v\ and 
V2, is given by: 

m 2 {v l7 V2) = [E{ Vl ) + E(v 2 )} 2 -[p{ Vl ) + p(v2)} 2 

= m 2 { Vl ) + m 2 (v 2 ) + 2[E(v 1 )E(v 2 ) - p(vi) ■ p(^)] ; (2.1) 

where bold quantities denote three-momenta. The maximum value taken by m{v\,V2) 
depends upon whether the a% particles can be brought to rest in the u rest frame. Such a 
configuration can be obtained in principle if the boosts of the a, particles in the respective 
8i rest frames can be arranged to exactly cancel the boosts they obtain in the to rest frame 
from the boosts of the 5$ rest frames. If this condition can be satisfied for some values of 
m{vi) then this configuration generates the maximum possible value of m(vi,V2)'- 

m max (v 1 ,V2) = m(u) - m(ai) - m(a 2 ). (2.2) 

The necessary and sufficient conditions imposed upon the masses by this condition are 
discussed below for the cases where the two 5$ decay chains are identical or different. 
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Figure 1: an phase space diagrams for toy models with identical decay chains. The hatched 
regions represent the annular rings referred to in the text. See text for details of models. 



2.1 Identical chains 

When 5\ and 82 are of equal mass, and similarly for a\ and ct2, the condition that both on 
particles be brought to rest in the oj rest frame is equivalent to the requirement: 



m Lx( w ) > m2 ( S ) ~ m(a)[m(u) - m(a)) > m mir », 



(2.3) 



where m max (v) and m m i n (v) are respectively the maximum and minimum possible invariant 
masses of all the visible decay products of each chain. If this condition can be satisfied 



then from Eqn. (2.2) 



m max (vi,v 2 ) = m(cj) - 2m(a). 



(2.4) 



In models satisfying Eqn. (2.3) the kinematic configuration saturating the bound Eqn. 
is that in which both a, particles are emitted against the motion of their parent <5j particles 
with sufficient momentum to cancel the boost in the oj rest frame provided by this motion. 

A useful means of visualising the kinematic configurations generated by a given model 
is provided by phase-space diagrams in which the cylindrical polar components of the 
velocities (3 of the two a, particles measured in the oj rest frame are plotted on the x axis 
(components parallel to the oj decay axis) and y axis (components perpendicular to the oj 
decay axis). Examples of such diagrams are shown in Figure [l]. In these diagrams the phase 
space accessible to each ai particle is represented by an annular ring, with the outer and 
inner boundaries defined by the kinematic configurations in which 771(1^) = ?r&min 

(vi) and 

m(vi) = m max (f) respectively. If m max (v) = m(5i) — m(a.i), as is the case for three-body 
5i decays, then the inner boundary possesses zero radius. This is because in this case the 
ai particles must be produced at rest in the 5i rest frames. 

Figure |l](left) shows the ai phase space diagram for a toy model satisfying Eqn. ([2.3|). 
The configuration saturating the bound Eqn. (2.3) is that in which (3 = for both ai 



particles, which requires that the origin lies within both annular regions. 

Figure |l](centre) and Figure [l|( right) show a» phase space diagrams for models which do 
not satisfy Eqn. ((2.3|), as can easily be seen by observing that in neither case does the origin 
lie within the annular regions. In such cases m max (vi,V2) is obtained from configurations 
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in which m(ai,a2) is minimised. These configurations must therefore simultaneously min- 
imise the velocity f3 of each a.i particle in the uj rest frame. In models similar to that shown 
in Figure |l]( centre) this requires that m{v\) = m(v2) = m m \ n (v), generating configurations 
lying on the outer boundaries of the annular regions. In Figure fright) this requires that 
m(vi) = m(v2) = m max (j)) generating configurations lying on the inner boundaries. In 
both cases the configurations saturating the bound generate /3 values lying on the f-axis. 
This can be understood by remembering that the closest point to the origin on a circle 
with centre displaced from the origin along a given axis lies at the intersection between 
that axis and the circle. 

In practice, if Eqn. ( |2,3| ) is not satisfied for a given model then m max (t>i, V2) can be 
determined from the maximum of the two values of m(v\,V2) obtained when m{v\) = 
m(v2) = rn m ^(v) and m(v\) = m{v2) = w m i n (f), assuming that both ai particles are 
emitted against the directions of motion of their parent 8{ particles. For fixed m{v\) = 
m ( v 2) = m(v) this configuration leads to the following expression for m(t>i, 7J2): 

m{vi,V2) = - — Trrpr m (0) - m (a) + m (v)+ 
2m z (d) L 

^— yj (m 2 (w) - 4m 2 (5))([m 2 (5) - m 2 (a) + m 2 (v)] 2 - 4m 2 (8)m 2 (v)) , (2.5) 



m{uj 



and so substituting respectively m(v) = m max (v) and m(v) = m m i n (v) into Eqn. (| 
enables m max (fi,t>2) to be determined. 

2.2 Non-identical chains 

When the two decay chains are not identical alternative kinematic configurations can max- 
imise m(yx,V2). In such models the phase space diagrams are not in general symmetric 
and qualitatively new configurations can be generated such as those appearing in the di- 
agrams in Figure In Figure §(left) and Figure ^(centre) the energy released in the 
decay of 5i is respectively less than or greater than that released in the decay of 82- In 
Figure ^(right) 171(82) is very much less than m(8\) leading to a larger boost of 82 and 
hence a larger offset of the centre of the left-hand annular region. It should be clear from 
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this diagram that the number of kinematic configurations which can in principle maximise 
m(vi,v 2 ) is greater than in the case of identical decay chains. 

As in the case of identical decay chains discussed above, the maximum possible value 
of m(vi,V2) is given by Eqn. (2.2) if the model can generate a configuration in which both 
on particles are at rest in the uj rest frame. The conditions for this to be the case for 
non-identical decay chains, analogous to Eqn. (|2.3[) for identical decay chains are: 



M > m 2 (5 1 ) - m( ai ) \m(uj) - m(ai) + "^"f^ 



^ j > m min ( Vl ), (2.6) 

To identify the possible configurations which can maximise m(vi,V2) when the condi- 



m 2 max {v 2 ) > m 2 (5 2 ) - m(a 2 ) m(w) - m(a 2 ) + ^^T^ 



tions Eqns. ( |2.6| ) and (2.7) are not satisfied let us first consider the energies and momenta 
of the oti and Vi particles in the Si rest frames. These may be determined from simple 
two-body kinematics to be: 



E( ai ) 
E{a 2 ) 



m 2 (5\) + m 2 (ai] 



2m{5 1 ) 
m 2 (5 2 ) + m 2 (a 2 ) 



m 2 (a\) + m 2 {v\) 



m 2 (v 2 ) 



2m(8 2 



E(v 2 ) 



m 2 (5 2 



2m{5 1 ) 

m 2 (a 2 ) + m 2 (v 2 ) 



2m{5 2 



(2.8) 



and 



p(vi) = p(ai) 
p(v 2 ) = p{a 2 ) 



v/^ 2 (ui) - m 2 (vi) = y / E 2 (a 1 ) - m 2 ( ai ), 



v 2 



m 



] (v 2 ) = yjE 2 (a 2 ) - m 2 (a 2 ), 



(2.9) 



where p(x) denotes the magnitude of the three-momentum of particle x. In addition the 
energies and momenta of the <5j particles in the uj rest frame are given by: 



E{5i 



m 2 {uj) + m 2 (5i) - m 2 (5 2 ) 
2m(w) 



E(5 2 



m 2 (uj) — m 2 (8i) + m 2 {5 2 ) 
2m(w) 



and 



p(5 l )=p(S 2 ) = ^E 2 {5 1 ) - m 2 (5 1 ) = y / E 2 (8 2 ) - m 2 {5 2 ). 



(2.10) 



(2.11) 



The magnitudes of the boosts ft provided to the <5j particles by the uj decay can be obtained 
from 

p(Si) 



ft 



E(5i 



(2.12) 



Now let us consider the possible extremal values of the momenta and energies of the 
aj and Vi particles in the uj rest frame, for fixed m(vi) and m(v 2 ). These can be obtained 
through Lorentz transformations of Eqns. (|2.8j ) and (|2.9| ) with ft values obtained from 
Eqns. (fnof)- (ED^) : 

E'(a\) = 7i[£ , («i) + c\Pxp(a\)\, p'(a\) = 7i[cif(ai) + ft-E(ai)], 
E'(a 2 ) = ^ 2 [E(a 2 ) - c 2 /3 2 p(a 2 )\, p'{a 2 ) = ^ 2 {c 2 p{a 2 ) - /3 2 E(a 2 )], 
E'(vi) = i\[E(vx) - c\fi\p{v\)], p'(yx) = 7i[-cip(«i) + fiiE(vi)], 
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E'(v 2 ) = j 2 [E(v 2 ) + c 2 fcp(v 2 )], p'(v 2 ) = j 2 [-c 2 p(v 2 ) - (3 2 E(v 2 )], (2.13) 



where primed quantities are measured in the uj rest frame and 7, is the Lorentz gamma 
factor derived from /3j. The quantities Cj = ±1 and are determined by the directions in 
which the a% particles are emitted in the bi rest frames. If we define b\ to be travelling in 
the +x direction in the uj rest frame then c\ = +1 corresponds to a configuration in which 
«i is emitted co- linearly with 5\, while if c\ = — 1 it is emitted contra-linear ly. If c 2 = +1 
then a 2 is emitted contra- linearly with 5 2 , while if c 2 = — 1 it is emitted co-linear ly. The Oj 
particles in these extremal configurations possess velocities given by p' (ai) / E 1 (on) . When 
m(vi) or m(v 2 ) are maximised or minimised these velocities define the coordinates at which 
the boundaries of the annular regions in Oj phase space diagrams such as Figure |] cross 
the x axes. The invariant masses of the visible particles in these extremal configurations 
can be obtained by substituting the energies and momenta of v\ and v 2 from Eqns. ( 2.13] ) 
into Eqn. ([O]). 

Now let us consider the kinematic configurations which generate extremal values of 
m(vi,v 2 ). We shall label these configurations with the signs of the Ci quantities, which 
denote the directions of the Oj particles in the 5i rest frames, together with the masses of 
the Vi particles. The notation we shall use is of the form {±„, ± w }, where the first(second) 
element inside the brackets refers to particle 1(2), the signs correspond to the signs of q, 
and the subscripts v and w describe the values of respectively m(v{) and m(v 2 ). So for 
example {+ m axi — mm} corresponds to a configuration in which a± is emitted co-linearly 
with 5\, a 2 is emitted co-linearly with 5 2 , m(v\) = m max (i>i) and m(v 2 ) = m m i n (v 2 ). 

Using our new notation, the possible configurations maximising m(vi,v 2 ) in the case 
of identical decay chains discussed in Section |2j] are {— m in, +min}, corresponding to Fig- 
ure |l](centre), and {- m ax,+max}i corresponding to Figure ||(right). In the case of non- 
identical decay chains we must also consider {— m ax,+min} and {— mim+max}- We have 
not exhausted the possibilities however. In order to maximise m(vi,v 2 ) we must min- 
imise both m(cei,a 2 ) and the net momentum of the ol\ol 2 system. In some models with 
highly asymmetric decay chains, for example that represented in Figure ||(right), these 
two requirements are mutually exclusive - decreasing 771(01,02) increases the net invisible 
momentum and vice versa. In such models we must find the values of m{v\) and m(v 2 ) 
which maximise m(vi, v 2 ) when respectively m(v 2 ) and m(vi) are maximised or minimised. 
These values are located at the turning points of m(vi,v 2 ) and are given by: 

m t P ( v i) = ™< 2 {8i) + m 2 (ai) 

m(a) 2m 2 (^)m(^)+^(^)+j3^) 

%) y/m^lm^Si) + m?{ Vj )} + m(5j)[AE{vj) + Bp( Vj )] ' 

where 

A = m 2 (uj) - m 2 (5i) - m 2 (5 2 ), B = ^ A 2 - 4m 2 (c5i)m 2 (5 2 ), (2.15) 



E(vj) and p(vj) are defined by Eqns. ( p^q ) and (|2.9|) , and j ^ i. Taking all possible 
combinations into account this leads finally to eight possible values for m max (vi,v 2 ) when 
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Eqns. fl2.6| ) and (2/7) are not satisfied. These values are obtained from the configurations: 



{ rniru "I'min}) { max) "I'max}) { 
{ tp) "I'min}) { tp) "I'max}) {" 



"miri) ~t~max}) { maxi "("min}) 
"mm) +tp}, {" "max) +tp}- 



(2.16) 



The value of m max _(vi, V2) is given by the maximum value of Eqn. ( |2.lD obtained by sub- 
stituting the energies and momenta from Eqns. ( 2.13| ) for the eight configurations listed in 
Eqn. (Ogl). 



3. Transverse mass end-points 

The transverse mass Mt of a set of visible and invisible decay products can be calculated 
from their transverse momenta and masses. The transverse momentum and mass of the 
aggregate visible decay product V of oj can be obtained by summing the four-momenta of 
the aggregate visible decay products x>\ and x>2 of each decay chain, while the transverse 
momentum of the invisible decay products is measured by the event E™ 1SS vector. The 
optimum definition of My depends upon the value of the lower limit x on the mass of the 
aggregate invisible decay product 1 . If this limit is zero, for instance because the m(aj) are 
unknown, then the optimum definition is: 

M |(0) = m\V) + 2[E T (V)Ef ss -p x (V) P r s -Py(V)pf ss }, (3-1) 

where 

E T (V) = ^p 2 T (V) + mi(V). (3.2) 

A configuration which maximises Mt(0) is that in which m(V) is minimised and pt(V) is 
maximised. This requirement implies that 81, 82, oi, 02, v\ and V2 must all move in the 
laboratory transverse plane. 

If the lower limit on the mass of the aggregate invisible decay product is non-zero, for 
instance because the m(«j) have been measured, then the optimum definition is: 

M|(x) = m\V) + x 2 + 2[E T (V)Et s \x)-p,{V) P t SS -Py{V)pf s % (3.3) 

where 

Ef ss { X ) = ^{E^Y + x 2 - (3.4) 

If the m(ai) are known then a conservative value for x is Tn(ai) + 771(02), which we shall 
use below. The absolute maximum value of Mr(m(ai) + 777(02)) is m(uj), which is obtained 
when the a% is at rest with respect to 02, and 8±, 82, ati, «2 5 v\ and ^2 are all moving in 
the laboratory transverse plane. Given the discussion of Section 0, we can use Oj phase 
space diagrams such as Figures [l] and Q to identify the configurations which satisfy this 
requirement. Note however that because we are now using only transverse momenta we 
must reinterpret such diagrams as representing the transverse cartesian components of 
velocity rather than cylindrical polar components about the oj decay axis. It is then clear 

Note that this is not directly a limit on the mass of the individual invisible decay products «i, for 
instance the LSPs in SUSY models. 
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that the configurations for which Mr(m(ai) + 771(02)) = 777(0;) are those located in the 
regions of the on phase space diagram in which the annular regions overlap, provided all 
the motion lies in the laboratory transverse plane. If the annular regions for a given model 
do not overlap, or some of the particles move out of the laboratory transverse plane, then 
this bound is not saturated. 



3.1 Identical chains 

In the case of identical decay chains Mt{0) defined by Eqn. fl3.1|) can be maximised with the 
configurations {+ m in> +min} and {— m in ; — min}> which generate in this case equal maxima. 
Substituting Eqns. ( 2.13j ) into Eqn. ( |3.1| ) and assuming that all the motion lies in the 



laboratory transverse plane we obtain the bound: 



M T ax (0) = m 2 (5) - m 2 (a) + ml lin (v)+ 



yj [m 2 (S) - m 2 {a) + m 2 min {v)] 2 - Am 2 (S)m 2 miri (v) , (3.5) 

where m m i n (7;) is the minimum value of the invariant mass of the individual aggregate 
visible decay products v\ and V2, as in Section |j| If 77i m i n (T>) = this reduces to 

M^{Q)=m{u)\l- T ^-\ (3.6) 

L 777^(o)J 

which is always less than 777(0;). 

In the case of Mr(777(ai)+m(a2)), i.e. M<r(2m(a)), the condition that M™ ax (2?77(a)) = 
777(0;) is equivalent to requiring 

777 2 (<5) — 777(a) [777(0;) — 777(a)] > 777^(7;), (3-7) 



which is a less stringent requirement than Eqn. (|2.3| ). Models satisfying this requirement 
include those represented in both Figure [l]( left) and Figure |I|(right). In the latter case the 
configurations saturating the bound possess en particles moving with (3 X = and /3 y / 0, 
in other words emitted transverse to the 5i direction in the uj rest frame. The fact that 
for a given model multiple configurations can saturate the bound rather than just one (as 
is the case for m(vi, V2)) shows that in principle the Mr(2777(a)) end-point can be more 
prominent, as shall be discussed in Section ||. 

If the requirement Eqn. ( |3.7l ) is not satisfied, as is the case for the model represented 
in Figure [l|(centre), then the configuration which maximises MT{2m{a)) is that which 
minimises 777(01,02), with all the motion in the laboratory transverse plane. This config- 
uration is {— m i n , +min}> which is also that configuration which maximises 777(7;!, ^2) given 



777(7;) = 777 m i n (7;) (see Section 2.1). In this case the aggregate invisible and visible transverse 



momenta are zero and thus from Eqn. 

Mf ax (2777(a)) = 2m(a) +m(v 1 ,v 2 ), (3.8) 
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or equivalently from Eqn. (2.5): 



Mf ax (2m(a)) = 2m(a) + ^$=r m\5) - m\a) + m' min (v)+ 



J (m 2 (w) - Am 2 (5)) ( [m 2 (5) - m 2 (a) + m^Jv)} 2 - 4m 2 (5)m^ in (f )) . (3.9) 



m(w) 

3.2 Non-identical chains 

In the case of non- identical decay chains the analysis for Mt(0) is very similar to that 
in the case of identical decay chains. Configurations which can generate the maximum 
values of M T (0) are again {+ m i n , ~l~min} and { — mim — mm}) although in this case they may 
in principle generate different values of M™ ax (0). Assuming that all the motion lies in the 
laboratory transverse plane, these configurations can be used with Eqns. ( gig ) and 



to generate two possible values of M™ ax (0) with the maximum of these two values used. 

When considering MT(m(cti) + m(ot2)) the analysis is more complicated. To simplify 
the discussion we shall first define decay chain 1 to be that chain which generates the 
smaller annular region in the on phase space diagram (see e.g. Figure ^). In other words 
we label decay chains such that the difference in f3 values for a± configurations + m i n and 
—min is less than the difference in (3 values for a<i configurations + m i n and — m i n . Given 
that we have already defined the +x direction to be that direction in which 5\ is emitted in 
the oj rest frame, this additional definition then maps Figure ^(centre) onto Figure ||(left). 

With this new definition of decay chains 1 and 2 we find that the requirement that the 
two annular regions in the on phase space diagram overlap is equivalent to requiring that 
the (5 value for the a\ configuration — m j n is less than the /3 value for the «2 configuration 
+ m i n , and the (3 value for the a% configuration + m i n is greater than the f3 value for the 
Oi2 configuration + ma x- 

These requirements are not satisfied by any of the three models 
represented in Figure |2| and in such cases the maximum value of Mr(m(ai) + 171(02)) is 
obtained from one of three configurations: {— mim +min} (as in the case of identical chains), 
{—max, - min} and {+mim +max}- The latter two configurations map to each other if the 
decay chains are defined as described above. 

4. Example: SUSY higgs decaying to gauginos at the LHC 

In order to demonstrate some of the simpler aspects of the above discussion a Monte Carlo 
study was conducted of heavy SUSY higgs particles decaying to the four lepton plus I?™ 1SS 
final state via X2X2 w ith three-body X2 lepton producing decays: 

H/A^x° 2 x° 2 ^UUxiXi- (4-1) 

The model chosen for study was the mSUGRA Point A model from Ref. Q, for which 
m(H/A) = 256 GeV, m(x§) = 110 GeV and m(xjj) = 60 GeV. The ISASUGRA 7.69 §1| 



RGE code coupled to HDECAY [22] was used to generate the sparticle mass spectra and 



branching ratios, while HERWIG 6.510 |23], |I| and ACERDET |2§] were used to generate 



and simulate 14 TeV LHC events corresponding to about 300 fb 1 of data (assuming 
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Figure 3: Distribution of detector-level four-lepton invariant mass values for H/A 
££££XiXi events from the Point A model. The dashed vertical line represents the expected end-point 
position from Eqn. (2.2). The small population of events lying beyond the expected end-point is 
generated by detector mis-measurement. 



perfect identification efficiency). Detailed studies of SM and SUSY backgrounds to the 
4£ + E™ ms channel analysis are beyond the scope of this paper - for this simple illustrative 
study only signal events were considered with a rudimentary detector level event selection 
requiring merely the presence of four isolated leptons in opposite-sign same-flavour pairs. 
The resulting values of m(v\,V2) = m(£i, £2, £4, £4), My(0) and My(2m(x5)) are shown in 
Figures || and |j. 

In the case of H/A — > X2X2 ~~ ^ ££££XiXi events the two 5{ decay chains are identical, 
with uj = H/A, Si = 62 = X2, ai = «2 = Xi an d i>i = «2 = £ + £~- In the example 
considered here m max (vi) = — m(xi) due to the three-body nature of the X2 decays, 

while m m i n (-Uj) = 0. Consequently the requirement in Eqn. ( |2.3D is satisfied and hence 
wimax(^i> ^2) = m(H/A) — 2m(xi) = 136 GeV. The requirement in Eqn. ( ^7| ) is also 
satisfied giving Mf? ax (2m(x?)) = 256 GeV. The M T (0) end-point is given by Eqn. ( p| ) 
leading to M™ ax (0) = 180 GeV. These expected end-point positions are represented in 
Figures || and || by vertical dashed lines and agree well with the observed end-points. It 
is interesting to note that the larger number of configurations saturating the bound on 
Mr(2m(x5)) compared with m(vi,V2) in this case leads to a more prominent end-point 
with a steeper gradient. This is true even at detector- level following smearing of the event 
E?p lss values used to calculate Mr(2m(x?))- 
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Figure 4: Distributions of detector- level transverse mass values for H/A 
events the Point A model. In the left(right)-hand figure the transverse mass is defined by Eqn. (3.1) 
(Eqn. (|3.3| )). The dashed vertical lines represent the expected end-point positions, given by respec- 
tively Eqn. and m{ui). The small populations of events lying beyond the expected end-point 
are generated by detector mis- measurement. 



5. Conclusions 

This brief paper has discussed the positions of end-points in the invariant mass and trans- 
verse mass distributions of the decay products of heavy particles decaying to pairs of 
semi- invisibly decaying products. The formulae presented here may prove useful for mass 
measurements if SUSY higgs bosons decaying to gauginos are observed at the LHC. The 
same techniques may also prove useful in other new physics scenarios, for instance given 
heavy states decaying via pairs of new W' bosons to massive stable right-handed neutrinos. 
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